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INTRODUCTION 


The research progress on calculations of interfacial stress in abrupt 
junctions and its effect on the defect state of the interface has yielded 
significant progress in the form of a series of papers. The first paper is 
coauthored with Professor J. H. van der Merwe and has been submitted to the 
Journal of Applied Physics. The second paper in the series is being complet- 
ed. As the title indicates the sharp interfaces in a multilayered super- 
lattice is treated for the case of equal elastic constants and layer thick- 
nesses. Part II of the series is for the same case but unequal elastic 
constants and thicknesses. Part III of the series is not well formed but is 
aimed at multilayers with different misfits such as those of step graded 
junctions or layers of altered lattice constants by means other than composi- 
tion variation e.g., variations in binding states from layer to layer. The 
series is expected to conclude with Part IV which will treat abrupt interfaces 
in mesa structures or other structures exhibiting finite width as opposed to 
infinite lateral extent. The following paper included in this progress report 
is the unreviewed part I of the series. 



AN EXACTLY SOLVABLE MODEL FOR CALCULATING CRITICAL MISFIT AND 
THICKNESS IN EPITAXIAL SUPERLATTICES. I. LAYERS OF EQUAL ELASTIC 
CONSTANTS AND THICKNESSES. 

Uan H. van der Merwe and W.A. Jesser 3 ^ 

Department of Physics, University of Pretoria, Pretoria 0002, 

Republic of South Africa. 

A parabolic interaction potential has been used to develop a model 
for calculating the misfit dislocation (MD) energy in the case of a 
superlattice of alternating layers of materials with equal elastic 
constants and thicknesses. The model, which is believed to be a good 
one for small misfits and to have some merit for covalent bonded 
materials, is exactly solvable for the critical thickness above which 
it is energetically favorable to lose coherency by the introduction of 
MD's into the interfaces. It was found, for a given misfit f, that the 
critical thickness for epitaxial superlattices free from their substrate 
is somewhat more than four times that for a single epilayer on a thick 
substrate. Furthermore, the critical thickness varies almost inversely 
with misfit to the power 1.22, when Poisson's ratio is 1/3. It was 
also shown that the critical misfit f obtained by equating maximal 
misfit strain and MD energies is a significant overestimate of f c „ 

The results for a superlattice are compared with those of a thin layer 
on a thick substrate. 

a) Permanent address: Department of Materials Science, University 

of Virginia, Charlottesville, Va 22901, U.S.A. 


2 


I. INTRODUCTION 

The concept of a limiting misfit as the upper limit of misfit below 

which epitaxial growth may occur, was introduced by Royer*. This 

2 

motivated Frank and van der Merwe to introduce the relevant concepts of 
misfit strain (MS), misfit dislocations (MD’s), critical misfit f and 
critical thickness h c> The critical misfit f £ is the limiting misfit 
below which a layer of fixed thickness is homogeneously strained into 
registry with the substrate and is the critical layer thickness 
above which an interface of given natural misfit loses registry by the 
introduction of MD’s. 

When epitaxy became involved in the fabrication of devices needing 

crystals of high perfection, knowledge regarding f £ and h £ assumed 

great importance and calculations have been made assuming equilibrated 

MD densities. Although notable successes had been achieved with this 

34 

approach * , observed densities were frequently greatly in excess of 

5-7 

the predicted densities . This discrepancy has often been attributed 

to non-equilibrium features of the growth process and to approximations 

needed to facilitate the analysis. Previous calculations have employed 

the Vo lterra approximation to MD's using the concept of a "cut-off radius" 

to account for the dislocation core energy, the presence of other MD’s 

5—8 

and the nearby free surface for films of finite thickness. While 

9 

the reliability of these approximations has been justified for the 

case in which the crystal thickness is of the order of the MD spacing 

and larger, this is not so otherwise. Furthermore an approximate 

. . 3 8 9 

procedure of energy minimization * ' has often been used too. This 
procedure, expressed by Eq. (42) is questionable particularly for small 
misfits. 
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The above regimes of concern, low misfit and thin crystals, are 
precisely those encountered in superlattice structures. It is thus of 
importance to resolve these issues so that confidence can be gained in 
the approximations adopted in calculations of f^ and h c » particularly 
in the case of superlattices . At the same time one can contribute 
towards solving the problems associated with the equilibrium accommodation 
of misfit by MD’s and MS for all misfits and thicknesses. 

In this paper an exact solution of the "thin film-small misfit" 

9 

case for superlattices is obtained using the socalled parabolic model . 

In this model the periodic interaction potential between the crystals 

across an interface is modelled by a succession of parabolic arcs rather 

than the sinusoidal description of the Frenkel-Kontorowa^ and Peierls 
1112 

Nabarro models * , Because of its sharply peaked crests it will be a 

more realistic representation for short ranged covalent bonding 
characteristic of many superlattices of semiconducting materials. 

In part I of this series we deal with the simplest case of a 
superlattice of similar materials where the layers have equal thicknesses 
as well as a thin film on a thick substrate of similar materials. This 
will be extended and generalized in subsequent parts to the more 
complicated cases of a thin superlattice film on a thick substrate and 
to superlattices of dissimilar materials and layers of unequal thicknesses. 

II. MODEL 

We make the conventional assumptions (i) that the crystals A (upper) 

and B (lower) have simple cubic structure with lattice parameters a and 
11 12 13 

b respectively ’ ’ , (ii) that a and b may be generated from a 

1 3 

reference lattice with lattice parameter c defined by 
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p = Pa = (P +l)b = (P + y)c, 


( 1 ) 


where P is an integer and (iii) that the crystals deform under applied 

forces like isotropic elastic continua with shear moduli y fl and y^ and 

11-13 

Poissons ratios v and V, . Equation (1) defines f, c and the HD 

Si D 

spacing p as 



ab 

c “ TTTT572 ’ 


_ c a - b 
p (a + b)/2 


( 2 ) 


The interfacial atoms are subjected to two competing forces: the 

atomic interaction across the interface which tends to align atoms on 

either side and the opposing elastic forces., The former is most simply 

modelled as a shear stress x having the periodicity c in the relative 

"displacement" U of interfacial atoms which were in registry in the 

11-12 

reference lattice . The two most common models for T are the 

11-13 

Peierls-Nabarro model that uses a sinusoidal dependence, and the 

9 

parabolic potential model in which there is a piecewise linear 
dependence between T and U: 


T (U) = y.U/d for |u| < d/2 = c/2, (3) 

where we may take the separation d of the two interfacial surfaces 
equal to c and the interfacial shear modulus y^ is a measure of the 
bonding across the interface. The relation (3) is repeated within 
every period c of U to form a saw tooth pattern. For mathematical 
convenience we consider at first a one-dimensional sequence of saw 
teeth as would obtain from one-dimensional misfit as defined in Eqs.(l) 
and (2). Two-dimensionality will be introduced in a subsequent section. 

The elastic relaxation introduced by t(U) generates an atomic 
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pattern with a sequence of localized disregis tries constituting a 

single sequence of misfit dislocations (MD's) * . They are spaced at 

intervals p and said to "accommodate" the misfit f. The misfit can 

also be partly accommodated by misfit strain e, i.e. a homogeneous 

strain which is superimposed on the oscillatory strains of the MD's, 

2-9 

to reduce or eliminate the misfit . Naturally this is only possible 
when one or both crystals are of finite thickness, which is the case 
in the crystal pairs constituting a superlattice. 

For mathematical simplicity the members of the pair will be assumed 
to have identical elastic properties. They will then be strained by 
equal amounts when they are equally thick. Thus, if the average lattice 
parameters in the strained configuration are a and b, the misfit strain 
e will be given by 


b -b a -a 
b a * 


( 4 ) 


It follows that in the strained configuration the relations (2) become 


(1 - e 2 ) - I - e 2 

P - P Z * c = c 


(1 - 2e/f ) 


1 - ef/2 


- = c = f - 2e 
p 1 - e 2 


(5) 


It is seen that f vanishes (p = °°) when e = -j f , as is implicit in 
Eqs. (4), and that to first order 


c = c and f = f + 2 e. 


( 6 ) 


The second relation in Eqs. (6) is the mathematical statement that the 
misfit is accommodated by MD's and misfit strain (MS) jointly. 

An important consequence of the presence of MD's is that opposing 
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interfacial surfaces may undergo different displacements w (x,0+) and 
(x,0-) normal to the interface where superscripts a and b relate to 
layers A and B respectively. Accordingly these displacements induce 
normal forces which we may, using Hooke 1 s law, model as an interfacial 
normal stress 

Ntw(x)i = ■ 2 -. i - [ w a (x,o +) r A x . >°~J l = .?. y ,L M ( 7 ) 

1 K }l l-2v. 1 d 1-2V. d • K ' 

i l 

where W(x) is the relative normal elastic displacement of the two 
opposing surfaces, \K an appropriate Poisson*s ratio of the interface 
and d = c. Figure 1 displays one period of a superlattice of alternating 
A and B layers, and the coordinate axes, with origin midway between MD's. 

We have introduced the simplifying assumption that A and B are 
elastically equivalent including the interface. 

This implies that 



= W, 


V. = V = V. = V 
i a b 


( 8 ) 


and allows us to have the layers equally thick; of thickness 2h * 2r)c, 
say. Apart from automatically satisfying the lateral force balance 
implied by Eq. (4), it also facilitates the simplifications in Eqs. (10) 
below. 

The bicrystal layer between midplanes PQ and P'Q* in Fig. 1 with 
the interface MN may now be taken as a representative unit of the 
superlattice. Careful consideration of the relevant physical principles 
and the geometry displayed in Fig. 1 guided us to adopt the following 
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simplifying assumptions for the interior bicrystal unit of an extensive 
superlattice containing MD’s: 

All field quantities (stresses p.j, T and N) have the periodicity 
of the MD spacing p (or p) e.g. 

P xx (x + p » z) = P xx (x,z) * (9a) 

Because of the source of antisymmetry, U(x), one also has antisymmetry 

in t(u(x)) of Eq. (3) and in p (x,0) = t(U(x)), hence p (x,z) is 

zx zx 

antisymmetric in x: 

p zx (-x, Z ) = - p zx (x,z). (9b) 

The stresses p and p , that act on the interface surfaces, are 
r zz r zx 

continuous across the interface, i„e. 


P zx (x,°+) = P zx (x,°-) -TlU(*)] 

P a (x,0+) = pNx.O-) -NlW(x)] . 
zz zz 


(10a) 


Because of symmetry p (x,z) and w(x,z) vanish on the midplanes PQ and 

ZX 

P’Q\ i.e. 

p (x,z) = 0, w(x,z) = 0 at z = ± h . (10b) 

Three more remarks are appropriate* (i) Because the interaction 
energy between the homogeneous misfit strain and the oscillatory 
strains that are associated with the MD f s vanishes they may be treated 
independently o This allows us, as in Eq. (9a), to use p for the MD 
spacing in the calculations and introduce the strained value p (see 
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Eq. (5)) when needed. (ii) Also, with the simple cubic crystals under 
consideration there will be a cross grid of MD's, rather than a single 
sequence. We make the simplifying approximation that the relevant 
energies of the HD’s are additive as is expressed in Eq. (31) below, 
(iii) The conditions (10) strictly apply to the interior of the 
superlattice, "far" (possibly two to four layers) from the free 
surfaces and "substrate" employed to fabricate the superlattice. 


Ill, GOVERNING EQUATION AND ITS SOLUTION 

We assume that the MD's are long and straight. The related problem 
is accordingly one of plane strain which is normally analysed in terms 
of an Airy stress function X(x,z) satisfying the biharmonic equation^ 


V"x = 0 

and, defining the stresses in terms of the relations 

n = i!x - - 3 2 x n= _a 2 x 

^XX ~ 2 ' ^ZZ ~ 2 * ^zx Sx3z 
3 z 3x 


(10 


02 ) 


The solutions of (11) satisfying the periodicity (9a) and having 

9 

the symmetry (9b) can be expressed in Fourier transforms 

X a ’ b = [ (A a ' b + zC a,b ) cosh mz + (B a ' b + zD a,b ) sinhmz]cos mx, (13a) 

m * 2-rrn/p, n = 1,2, 00 (13b) 

where A, B, C and D are as yet unknown Fourier coefficients and 
summation is implied in all analogous subsequent expressions. 

By applying (12) we obtain results of the form (deleting superscripts): 
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{ x* } l"c) +2D,Coshmz * 2C )sinh»zUosmt, (14) 

z z 


P zx ~ m f (®(B +zD) +C}cosh mz+ {m(A +zC) +D}sinh mz] sin mx. 


In order to employ Eqs. (10) we need expressions for the x- and z- 
components of displacement u and w. We first obtain from Hooke's law 
the strains e.. as 


2we xx ■ 1 (l - v) Pxx - ■ ue 2 x ■ p zx : 


05 ) 


7uo _ m r / m(A + zC) +2(l-v)D, , . m(B + zD) +2( l-v)C, . , , 

2lJe { xx } m -m(A + zC) - 2VD > COsh raZ+ { - m (B +zD) -2vC }sinh mzlcosmx » 


U e zx " (ra(B + zD) +C}cosh mz+{ra(A +zC)+ D}sinh mz] 


sm mx. 


06 ) 


We now calculate consecutively the rotation and the displacements u 
and w by path integration of the form 


u " u o + /( li dx * 7? dz) ' u o ■ °' 


where 


9w 

—i l = 

9x 



9e 9w 9e 

xx 2 . _ 1 zx 

9z * 9z ~ 2 ~Tz 



3vi _ 3u _ 1 . .. 9w _ 1 9w 

9x xx' 9z 2 zx y* 9x 2 e zx W y’ 9z " e zz 


07 ) 


08 ) 



10 


We obtain 

2lJWy = 2(l-v)m 2 [D coshmz +C sinhmz]sin mx . 

2yu = [ {m(A + zC) +2( l-v)D}cosh mz+{m(B +zD) +2( l-v)C}sinh mz] sinmx (19) 

2pw =~[ {m(B + zD) ~ ( l-2v)C}cosh mz +{m(A + zC) - ( 1— >?v)D}sinh mz] cos mx 


and hence, at z =0 


v xd . a b xd L r m / .a .b\ 1-v , a b. , . 

U(x) = — + u - u = — + 1 ^ (A - A )+ — — (D - D ) ] sin mx, 


1 




W(x) = ~ [mCB 1 ' -B a ) + (l-2v)(C a -C b )] cosmx, 


( 20 ) 


where we had used Eqs, (8), and xd/p is the nonelastic vernier relative 
displacement due to the misfit. 

The governing equations, i.e, the equations defining the Fourier 
coefficients are obtained by substituting from Eqs. (3), (7), (14), (19) and 


(20) into (10): 

-m 2 A a + m 2 A b = 0 (a) 

m(mB a + C a )-m(mB b + C b ) * 0 (b) 

m{m(B a + hD 3 ) + C a }c + m{m(A a + hC 3 ) +D a }s = 0 (c) 

m{m(B b - hD b ) + C b }c -m{m(A b - hC b ) +D b }s = 0 (d) 

(21) 

{m(B a + hD a )-(l-2v)C a }c +{m(A a + hC a )-( l-2v)D a } s = 0 (e) 

{m(B b - hD b ) - (l-2v)C b }c -{m(A b - hC b )-( l-2v)D b }s = 0 (f) 

-m 2 A a + -^jyd {mBb - (1 " 2v ^ b} = 0 (g) 

, B a a. m , .a A b. 1-V , n a _b. y(-l) n+I ... 

m(mB + C ) - (A - A)- -j- (D - D ) = — , (h) 



J] 


where 

c = cosh mh and s = sinh mh , (22) 

we had used the elastic equivalence of crystals A and B (See Eqs. (8)) 
and (-l) n+ *d/irn for the Fourier sine coefficient of xd/p. Note that c 
in (21) and(22) is not the same as the reference lattice constant c in 
Eq. (1). The latter we also approximate by d. This system of linear 
equations may be solved to give 


A* 5 = A 3 = 0, A = -m 4 (mh + sc) - 2ym 3 s 2 

B b . B a _ _ _ D a D a _ y(-l) n+1 

B -B 2 " ~~2 ’ ~K VnlP' +Q7 

c c 

(23) 

P = (mh + sc)m 4 c 2 [m(mh +sc) +2ys 2 ] 

Q = 2(l-v)d 1 m 3 [m(mh +sc) + 2ys 2 ] . 
where 

hy = 2(l-v)h /(l-2v)d - An for v = 1/3, n = h/d = h/c . (2A) 

Important quantities which we need for subsequent calculations are 

p (x,0)and U(x). We obtain by substituting from Eqs. (23) into (1A) 
zx 

and (20) 

oo 

p zx (x,0) = 2 P^D a s * nnut » (25) 


U(x) = — + 2 Q D smmx, 

p . n n ’ 

K n=l 


( 26 ) 


where P > 0 and D 11 are the P, Q and D 3 in Eqs. (23). 
n n n 

Also the quantity Q/P plays a significant role and may be 
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written as 


*<« 


cs + X 


where 
8 = 


^ _ f J!_\ r 

TT^vTp " W f 


and 


X = 2ithn/p = 2irnfn = £n 


( 27 ) 


(28) 


and Yh is defined in Eq. (24). 

The relations (23) - (28) now constitute the solution of the 

problem concerned with the atomic relaxations resulting from the 

competing interactions. A plot of U v.s, x will display the resolution 

9 13 

of the interface into a sequence of MD's spaced at distances p * 

When MS coexists with MD's p,c and f must be replaced everywhere by 
the MS dependent quantities p, c and f defined in Eqs„ (5). 


IV. ENERGY 

The energy per unit area of an interface consists of several 

contributions: the MD strain energies E^ and E^ in the two crystals and 

the energy of misfit E due to the residual disregistry at the interface. 

rn 

The homogeneous misfit strain energies E_ and E in the two crystals may 

5 e 

simply be added to the above energies to yield the total energy. The sum 

3 15 

of the latter two energies is given by the well known expression * 


E- = 4[ (l+v)/(l-v)]ycne 2 (29) 

e 

which incorporates the equality of the elastic constants (Eq. (8)) and 

equality of MS along the interfacial cubic axes. It has also previously 
been 
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shown that the energy per unit area of interface associated with the 
MD ' s i*e. Ep = + E^ + E^ in units of yc/4iT 2 , is given by 


= e d / 


_ J_ f p/2 
2 2p J -p/2 


j cx 
dx — 


4tt 


P 2x (x, 0) / 


yc 

2 


4tt 


00 

2 


1 


n=l nln + 8 ^(X)] 


(30) 


where the second line follows on using Eqs. (25) - (28). ' Since in the 
case of quadratic symmetry under consideration there is a cross-grid of 
MD’s the total energy per unit area of the interface is given by 


E - E 5 * 2 h - - - ‘ ■ w cn 1 


e 

7 


1— v 


16 m 2 (i+v)n 


s D (M)I, 


(31) 


where everywhere 8 and X should be replaced by 8 and X i.e. f by f, to 
stress their dependence on the MS e as defined by Eqs. (5) and (6). 

We can execute two consecutive approximations to the results in 
Eqs. (30) and (31), which are of interest. The first one is the thick 
(h,n -*■ «*) limit Ep. It follows from (27) and (28) that in this case 
(provided f f t 4 0) . 


4>(X) = 1 for n *+ 00 , and 


oo 


1 


°° _ CO JJC_ _ y 

D " D 47r 2 n=1 n(n + 


(32) 

(33) 


• • • J 

This agrees with a previous result and may be written in terms of 
derivatives of the gamma function T(x) 16 as 


e* = 6U(1 +1/6) + Y I , 


(34) 
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where ip(x) is the digamma function r T (x)/T(x) and y ** 0.5772 .... is 
Euler's constant. 

The second approximation to the result in Eq. (33) is to replace 
the sum by an integral: 



dn 

n(n + 1/g) 


$ Hn(l + 



(35a) 


Both the value of the lower integration limit n^ and the accuracy of 
this integral in representing the sum may be assessed through a 
comparison of the relations in (34) and (35). As n = 1 for the first 
term in the sum of Eq, (33), the appropriate lower limit of the 
integral in Eq. (35a) is expected to be somewhat less than unity. 
We obtain the value of as 


= exp(-^) = 0.56! 


(35b) 


by comparing Eqs. (34) and (35a) for large values of 1/0. 
In this case 


Ep - 6 InieYjV) «* 6 £n(l/n.0), 


(36) 


which gives the result in Eq. (35b) for n.. 

l 

We now return to the general case of Eq. (30) and adopt the 
approach employed for the "infinitely thick" layer in replacing the 
sum by an integral: 


— CD 

C/ V 


dX 

x[x + a<Kx)l 


CI(n.£,a). 


(37) 
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Using the transformations and definitions of Eqs. (28) 
we obtain 


X = mh = £n; Z, = 2nT|f, n - h/c , f = f - 2e, 

(38) 

c/B = 2(i-v)n = ot , 

where the "bar" signals the presence of MS e. 

Near the critical thickness, which is presently an important 
regime X 0 because f 0 and hence <f>(X) °°. The series accordingly 

converges exceedingly slowly and the integral approximation (37) becomes 
more convenient. However, the integrand is rather complicated and does 
not lend itself to exact integration. Since <j>(X) tends to unity as X 
becomes large and is already almost unity at X * 2 it remains to find 
an acceptable approximation for small X. Power series expansion of 
c (= coshX) and s suggests a sufficiently simple representation of the 
form +<X) = A/X + BX where A and B are rational functions of f|. 

Numerical comparison has shown that 


<J>(X) = 


J 2X 


(1 +^) 
1 


for X <2 
for X >2 


(39a) 


V 

is an acceptable approximation. It is accordingly convenient to 
divide the interval (n^C» °°) of integration into two subintervals 
(n^C, X t ) and (X fc , 00 ) where the transition value of X is 


X«. = 2 


+ n-C 

i 


(39b) 


We obtain on integration 
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e D = 5I(n-0 = Ji{-| [ arctanBX t - arctan Bn^cl + £n(l + JL)}, (40a) 

X t 


where 


;,Bc 

on 


(40b) 


and the n.£ dependence of I(n.£) relates only to the lower limit n.E 
1 1 1 

in (37) and not n.£ in (39b). 

A comparison of the integrated relations (40) with the sum in 
Eq. (30) (summing over sufficient terras to reach (n = X = 2 and 
integrating from X * 2 onwards) shows that the discrepancy does not 
exceed 6%, The sum is the greater of the two. Note, however that 
the relations (40) should be restricted to misfit values f < 1 /ttt). 

Above this value Eq. (35a) should be used. 

In the limit £ = 27Tr|f 0 (f -*■ 0) the thickness approaches the 
critical value and iCn^C) in (40) becomes 

K°) - £§■ arctan 2B + i £n( 1 + ^ ) . (41) 


When however, 0 and hence £ = 2irnf 2TiT|f becomes large at nonzero 
f, e D as defined by (40a) approaches the value in Eq. (35), as it 
should. 


V. EQUILIBRIUM: CRITICAL MISFIT AND THICKNESS 

In the past the critical misfit f at which the minimum energy 
configuration undergoes a transition from a coherent interface to one 
with MD's has often been approximated by setting the energy for 
the coherent configuration (MS e = j f for a superlattice) equal 
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to one with only MD's, characterized by natural misfit f c , i.e. 

V f c’ n) = 2 V f c’ n) ’ (42) 

the contention being that there are either no MD's or enough to 

accommodate practically all the misfit as in the one-dimensional Frenkel- 
2 

Kontorowa model . Alternatively, the critical thickness h = n c is 

c c 

the thickness at which equilibrium considerations require MD's to be 

injected for given misfit f. Thus n is defined by (42) by setting 

c 

f = f and n = t) . 
c c 

The true criterion for the minimizing strain e^ is 

3E/9e =0 at e = e . (43) 

m 


Among others we wish to assess the accuracy of (42) „ Since the integral 
formulation of e D in Eq. (37) facilitates the numerical evaluation of 
Ejj we perform the minimization of E in Eq. (31) using the relation (37) 
rather than the series in (30). We obtain after some simplification 


[I( V> *«. * 

where 


X 

m 


n.C , 

l m 


m 


27171? = 27in(f - 2e ) , 

m m 


(44a) 


(44b) 


as may be seen with reference to Eqs. (6) and (38). Hence, for a given 
thickness ru the critical misfit f follows by letting e f /2 
yielding 


( 45 ) 
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interface is large, one needs analytical expressions that accurately 
give the energy of the MD arrays per unit area per interface. Both 
of these regimes of interest have been satisfied through the 
development of Eqs. (37) - (41). It is instructive to compare the 
energy values obtained from this result for different misfit values, 
which must be small because of the approximations involved in Eq. (6). 

Figure 2 displays the functional dependence of the MD energy 
per unit area per interface of a superlattice free from its substrate. 
The thick curve labelled °° is a plot of Eq. (35) which represents 
"infinitely thick" layers. The numbered curves are plots of Eq. (40) 
with each number representing the thickness value. For purposes of 
comparison the MD energy for a thin layer on a thick substrate is also 
shown. The superlattice curve labelled 1 closely approximates 
the curve for a monolayer on a thick substrate up to values of f <0.1. 
The thick curve labelled °° also represents the asymptotic value of 
the thin layer on an infinitely thick substrate for the layer 
thickness being large. Because of the choice of a representative 
section of the superlattice for purposes of calculation, a monolayered 
superlattice corresponds to n = (see Fig. 1) and is labelled by 
this thickness value. Only for the infinite layer is the slope of 
the curve infinite at the origin. In all cases the slopes detrease 
with increasing misfit at constant thickness and also decrease at 
zero misfit as the thickness value is decreased. The nearly linear 
rise of the energy with misfit for a very thin layered superlattice 
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can be understood from the arrays of MD's which act nearly 
independently of their in-plane spacing and therefore contribute to 
the superlattice energy essentially in proportion to their number 
density. 

One notices that there are regions of f for which the value of 

OO 

Cp exceeds that value for the infinitely thick layered 

superlattice. The misfit values for which this occurs are near 
nf = 1 /it which is also near the transition value of Eq. (39). 

OO 

The rise of over is not caused by the approximations 

involved in Eq. (39) but is also found in the exact representation 
Eq. (30). Mathematically it occurs because <f>(X) in Eq. (30) falls 
below unity near X fc while $(X)* unity for all X when n “*■ 00 
(Eq. (33)). Physically this result can be understood as a consequence 
of the boundary conditions for the superlattice case. There is the 
constraint of zero displacement at the midplanes of a superlattice 
layer while the free surface of the thin layer on a thick substrate 
is stress free. 

A. comparison of the MD energy of a superlattice with that of a 
single layer on a thick substrate shows that when the layer 
thickness of the superlattice is small, strong cancellation of the 
stress fields of the MD's in adjacent layers occurs provided the 
misfit is not too large. The magnitude of this cancellation 
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reduces the MD energy by about 20% in a tnonolayered superlattice of 
Small misfit as compared to a single monolayer on a thick 
substrate. 

Also shown in Fig, 2 is a graphical comparison of the two methods for 

determining at a fixed thickness q, the critical misfit f above which 

c 

it is energetically favorable to introduce MD's into the interface of a 
superlattice. This pair of misfit-thickness values is exactly equivalent 
to giving the critical thickness n c (-n) for loss of coherency in a 
superlattice of fixed misfit f(=f c ). The one method of equating the 
homogeneous strain (MS) energy of a coherent superlattice to the MD 
energy when all the misfit is accommodated by MD's is given analytically 
by Eq. (42) and graphically by the intersection point of the homogeneous 
strain energy (plotted against f ■ 2e) with the MD energy curve. Plots 
of this situation are shown for two thickness values. The homogeneous 
strain energy curves A and B (dotted parabolas) correspond respectively 
to half-layer thickness values of q equal to 10 and 1. The intersection 
points are near misfit values respectively of 0.0315 and 0.154 which are 
also the calculated values from Eq. (42). The second and proper method 
of calculating the critical misf it-(critical) thickness pair of values 
is equivalent to equating the slope of the MD energy curve at the 
origin (all of the misfit is accommodated by MS) to the slope of the 
same value on the homogeneous strain energy curve. This is the 
condition expressed by Eq. (43) for = -I f a nd is shown in Fig. 2 for 
the case q = 1 by the two (dashed) tangents T, One tangent touches the 
MD energy curve at the origin to determine its slope, and the other 
tangent of the same slope is placed in contact with the homogeneous 
strain energy curve at the unique point of the same slope labelled i. 
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This point is at the misfit value 0.131 which is also the calculated 

value from Eq. (45). Notice that the misfit values given by the two 

methods, one proper and labelled i, the other approximate and 

labelled ii are rather different (0.131 versus 0.154); ii is too high 

by about 17%. Clearly the proper method is preferred and offers an 

. _ 3,8,9,17 

improvement over the approximate one. 

An exploration of how the magnitude of this overestimate varies 
with thickness revealed that it varied from about 17% for a mono layered 
superlattice to about 35% for large values of n, i.e. over a thousand 
monolayers. A similar comparison of the estimated critical thickness 
for a given misfit showed that the estimate exceeded the proper value 
by about 30% when the layer thickness is small and by about 40% when 
ri is large. This result suggests that one could "correct" the 
estimated value of critical thickness by multiplying it with a factor 
of about 0.75. 

The main objective of the present considerations is to find an 
exact expression for calculating the critical thickness n c above 
which loss of coherency at a given misfit is energetically favorable 
(Eq. (45)). The evaluation of this expression for the case of uniform 
elastic constants (y g = y^ = y^ = y) and v = is displayed in the 
curves of Fig. 3 0 Here one sees a strong dependence of n c on 
natural misfit f for small values of f. From a log-log plot (inset 
on Fig. 3) of the thickness-misfit values a nearly straight line 
emerges so that over a range of values a power law models the 
dependence fairly accurately. It follows from the figure that 
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n c * 0.114 f' 1,22 , (46a) 

or equivalently 

f » 0.17 n'°‘ 82 » (46b) 

These expressions describe the equilibrium case for which there are 
no barriers to the acquisition of MD’s. In practice barriers exist 
and therefore one expects significant departures from these 
equilibrium values when comparisons with experimental values are made. 
These differences between experiment and theory are common with the 
experimental values of MD densities exceeding the theoretical ones 
particularly in the case of semiconductor materials where the frictional 
force opposing the motion and nucleation of MD’s can be quite large. One 

can see the magnitude of these differences for example in the work of 

3 8 

Matthews and People and Bean . 

It is of interest to compare the results of the present development 
for superlattices of individual layer thickness 2r\ to the same parabolic- 
potential model but for an infinite substrate and a single epitaxial 
film of equivalent thickness 2n (Appendix)* There are two differences 
between these two cases* The first difference is that for the case of 
superlattices, the parallel MD arrays in adjacent interfaces strongly 
interact when the layer is thin and lower the energy of the MD array 
in comparison to that of the single-interface case on the infinite 
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substrate. The second and more important difference lies in the 
accommodation of misfit by homogeneous elastic strain. For the 
superlattice each layer accommodates half the misfit so that the N 
interfaces and N+l layers can become coherent with one half of the 
elastic strain that is required to produce a coherent interface in the 
infinite-substrate case. Because of the square dependence of the 
homogeneous strain energy on strain, this introduces a factor of four 
between the two cases. The combined effect on critical thickness of 
these two differences is displayed in the Table. 

The thickness of the single layer is compared to twice the half- 
thickness of each superlattice layer. The ratio of layer thickness for 
loss of coherency in the case of a super lattice, as compared to that of 
a single layer on a thick substrate is about five for small layer 
thicknesses and approaches four for large thicknesses. This result 

agrees with factor of four argued for superlattices by Matthews and 
1 8 

Blakeslee . The ratio of critical misfits for these two cases also 
approaches four for thick layers but starts at a factor of about three 
for a thin layer. 

The ratio of critical thicknesses greater than the factor four 
introduced by the homogeneous elastic strain term can be understood in 
terms of the slope of the MD energy as a function of misfit. The MD 
energy for the superlattice case is less than that for the case of an 
infinite substrate. Correspondingly the slope of this curve is also 
less at the common value of misfit at the origin. In satisfying the 
parallel tangent condition this results in an overall shift in critical 
thickness values by a factor somewhat more than four when comparing 
the superlattice case to the case of an infinitely thick substrate. 
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VII. CONCLUSIONS 

A parabolic interaction potential has been applied to the case of 
a free standing superlattice of many layers to obtain an exact solution 
for the energy of the arrays of misfit dislocations at the interfaces. 

An analytic expression for this energy as a function of misfit and 
layer thickness is provided in Eqs. (37) and (40). The main results of 
this theoretical development are: 

1 . The energetically favored critical thickness for loss of coherency 
depends approximately inversely on misfit to the power 1.22 for the 
case of uniform elastic constants and a Poisson's ratio of 1/3. 

2. A comparison of the critical thickness in a superlattice with that 
of a single layer on a thick substrate shows that for a fixed misfit, 
the critical thickness of the superlattice is between four and five 
times that of the single layer. This result occurs because for the 
superlattice case each layer accommodates half the misfit and the 

derivative 3e D /3f of the MD energy, is somewhat smaller than that of 
the thick-substrate case at f 0, except for large layer 
thicknesses where the values of approach one another. 

3. Estimating the critical thickness and critical misfit by equating 
the maximum strain energy (e = f/2) to the maximum misfit dislocation 
energy (e = 0) was shown to be rather crude; the estimated critical 
misfit for given thickness is too high by about 25% for all 
thicknesses. Similarly the estimated critical thickness for given 
misfit is too high by about 35%. 

4. The equilibrium elastic strain and hence the lattice constants of 
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the superlattice vary nearly inversely with thickness once the 
critical thickness for loss of coherency has been exceeded. 

5. The pointed high crests of the parabolic potential, which is the 
main objection to the model has on the one hand little impact when 
the misfit is small (MD’s far apart) and otherwise may have merit 
for short ranged covalently bonded materials. Its main advantage 
is that it yields an exactly solvable model for the critical misfit 
and critical thickness. 

6. The MD energy per unit area per interface for a thin layered 
superlattice can rise slightly above that for an infinitely thick 
layered super lattice when rjf is near 1/tt. 

7. The analysis applies to the interior of the superlattice far from 
the free surfaces. 
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APPENDIX: THIN EPILAYER ON A THICK SUBSTRATE 

A comparison of the thin epilayer -thick-substrate case with the 
superlattice had always been of great interest. This interest still 
persists. The former can be dealt with by some modification of the 
superlattice analysis. Although this had been dealt with previously 
the results needed for the present comparison have not been given and 
can be derived rather compactly from the foregoing calculations. 

Let A be the thin layer of thickness h and B by a semiinfinite 
substrate. The main modification in the model is that the midplane 
boundary conditions (10b) need be replaced by 

p ® x (x,z) - 0, p* z (x,z) = 0 at z = h (A. 1) 

for the free boundary of A and 

p..(x,z) ->0 as (A. 2) 

for the vanishing of the stresses in B far from the interface; i (and j) 
being either x or z. 

The condition (A. 2) can be incorporated by deleting the e mZ parts 
of X \ (Eq. (13a)) and in subsequent expressions relating to B. This 
is most simply handled by putting 

A b = B b , C b = D b . (A. 3) 


We may thus write down the quantities 
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b r /*b . ~b N , h, mz . 

p =mlm(A + zC ) + C ] e sintnx 

x z 

b 2 / »b -b* inz 

p = - m (A + zC ) e cosmx 
zz 

(A. 4) 

u b (x,z) =-j|j[m(A b + zC b ) + 2(l-v)c b ] e mZ sinmx 
w b (x,z) = - -^|j(m(A b + zC b ) - ( l-2v)C b ] e mz cosmx 

which we need explicitly below and in which we had used the elastic 
equivalence in Eqs. (7). It now follows from Eqs, (19) and (A. 4) that 

U(x) = xd/p + {mA a + 2(l-v)D a - mA b + 2(l-V)C b }/2y 

„ a Vi u (A « 5) 

W(x) = {-mB + (l-2v)C + mA -(l-2v)C b }/y . 


Substitution from Eqs. ( 14) , (A.4) and (A. 5) into (10a) and from (A. 4) 

into (A .1) shows that Eqs. (21) (a) - (c) and (h) are unchanged and 

that otherwise we obtain for p a (x,h) = 0 in (A.l) and p a (x,0+) = N[W(x)] 

z z z z 

in (10a) the relations 


-m(A a + hC a )c - m(B a + hD a )s = 0 

-m 2 A 3 + - 77 -y [ mB a - (l-2v)c a - mB b +(l-2v)C b ] = 0. 

\ l-zV;d 


(A. 6) and (21) (a) - (c) and (h) are six equations needed to solve for 
the six unknown Fourier coefficients. Our main interest is to 
calculate using the integral in Eq. (30). Hence we only need the 
values of B 3 and c in 


zx 


m[ mB 3 + C 3 ] sinmx. 


(A. 7) 


We obtain 
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B 3 = -(mV + <5mV)R/A, 

C 3 = [ m 2 s 2 + <5m(sc + s 2 + mh)]R/A f 
A = m(P + Q) , R = (y/irn)(-l) n+1 , 
p = m 3 (s 2 - m 2 h 2 ) + m 2 <5(s 2 + sc + mh - m 2 h 2 ), 

Q = Xm 2 (s 2 + cs - mh - m 2 h 2 ) + mX6 (c + s) 2 , 

X = (l-v)/d, 6 = X/(l-2v), a •III t 

p (x,0) = R/(l + QP" 1 ), 

4A 

00 __ 


P = 2 1 — „oo A Y 

D ! I = C / r — 

ntn + <t>(x)/B] P X[X + a<|>(X)] 
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♦ <x) = 


mQ = BnhQ = \_ X(s 2 + cs -X - X 2 )-*- 6h(c +s) 2 
2XP hP 2 X(i 2 -X 2 ) + 6h(s 2 +sc +X-X 2 ) 


9 


(A.8) 


where c, s, and X have been defined in Eqs. (22) and (28). At this 
point it is of interest to note that 4>(X) 4 I as X 4 ” (r) 00 ) as 

before so that the same result as in Eq. (33) is obtained for large 
thicknesses but that no simple expression like (39) for small X, could 
be obtained that connects accurately enough to <t>(X) = 1. 

With respect to energy minimization we obtain, instead of the 
relations (6), (31), (38), (44a) and (45): 


f C = f t + e t , = 2im t (f t -e), X C =nC fc , X fc -n.H , a* =2(l-v)n t 

m l Til 


e . ( ' . - v > ■ c t i(ii t ,o t 

1 "V ^ 0^2/1 ^ 


)] 


8tt (i +v)n 


(A. 9) 




with I(X t ,a t ) and 1(0, a*") given by Eq. (37) not by Eqs. (40a) and (41), 
m 
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where the letter t designates (t)hin f ilm-(t)hick substrate. It is 

notable that the expression for e fc in (A. 9) is of the same form as e 

ra tn 

in Eq. (44a) but that the expression for f^ in (A. 9) differs from 

that for f in Eq. (45) by a factor of two; that is apart from the 
c 

fact that <}) t = <J> in Eq. (A. 8) is a different function and that we 
usually make a comparison where the layers in the superlattice have the 
same thickness as the thin film, i.e. x\ 
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TABLE, Comparison of the critical thickness v\^ of 

c 

a single layer on a thick substrate and 2r\ of a 

c 

layer in a superlattice when the natural misfit f is 

the same in both; similarly the critical misfits f 

and f* respectively are compared for equivalent 

thicknesses ri** 52 2r| • 
c c 


f fc 

H 

2 \ 

2l c /tl 

.07143 

l 

5.09 

5.09 

.04395 

2 

10.03 

5.04 

.01283 

10 

49.3 

4.93 

.001955 

100 

460 

4.60 


t 

c 


V £ 


t 

c 


2.78 

2.99 

3,45 

3.66 
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Figure Captions 


Fig. 1 • 


Fig. 2. 


Schematic representation of a superlattice of alternating 
layers A and B of equal thickness 2h = 2nc. The AB and BA 
interfaces are represented by solid lines and contain 
planar arrays of MD f s of alternating sign as shown. The 
coordinate origin is in an interface and lies midway between 
MD f s with the z-axis normal to the interface planes and the 
y-axis along the axis of the MD's shown. A representative 
section of the superlattice is given by the mid planes r f Q ? 


in layer A and PQ in layer B each lying at a distance h = r|c 
on either side of the reference interface MN that contains 
the origin and is located in the interior of the stack of 
layers . 

Misfit dislocation energy per unit area (numbered solid 

curves in units of Pc/4it 2 ) as a function of misfit f and 

homogeneous elastic strain energy per unit area E_/ (yc/ 2tt 2 ) , 

(lettered dotted curves) as a function of e = f/2. For 

calculational purposes y = y - y. = y and v = 1/3. The 

a d i 

thick curve labelled 00 corresponds to Eq. (35) while the 

rest of the numbered curves are given by equation (40) with 

the numbers referring to the half- thickness r| = h/c. 

Strain-energy curves A and B refer respectively to n = 10 and 

1. The two linear tangents T are drawn by dashes for the 

case n = 1 and define the critical misfit f c by 3(2E^+ E-) /3e =0 

at 5 = f /2 given by point i. Point ii represents an 

me 

approximation to the critical misfit as calculated by 

2E n (f ) = E~(e = f / 2). Note that the value of critical 
u c e m c 
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misfit at i (f ' 0.131) is rather different from that of 
c 

point ii (f c ~ 0.154). 

Fig. 3. Graph of the relationship between critical thickness 0 C and 

the corresponding critical misfit f c> Note from the inserted 
log-log plot that the critical thickness is reasonably 
approximated by a power law with an exponent near 1.22. 
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